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ABSTRACT 

Simulations of the formation of Population III (Pop III) stars suggest that they were much more 
massive than the Pop II and Pop I stars observed today. This is due to the collapse dynamics of metal- 
free gas, which is regulated by the radiative cooling of molecular hydrogen. We study how the collapse 
of gas clouds is altered by the addition of metals to the star- forming environment by performing a series 
of simulations of pre-enriched star formation at various metallicities. To make a clean comparison to 
metal-free star formation, we use initial conditions identical to a Pop III star formation simulation, 
with low ionization and no external radiation other than the cosmic microwave background (CMB). 
For metallicities below the critical metallicity, Z cr , collapse proceeds similarly to the metal-free case, 
and only massive objects form. For metallicities well above Z cr , efficient cooling rapidly lowers the 
gas temperature to the temperature of the CMB. The gas is unable to radiatively cool below the CMB 
temperature, and becomes thermally stable. For high metallicities, Z > 10~ 2 - 5 Z Q , this occurs early in 
the evolution of the gas cloud, when the density is still relatively low. The resulting cloud-cores show 
little or no fragmentation, and would most likely form massive stars. If the metallicity is not vastly 
above Z cr , the cloud cools efficiently but does not reach the CMB temperature, and fragmentation 
into multiple objects occurs. We conclude that there were three distinct modes of star formation at 
high redshift (z > 4): a 'primordial' mode, producing massive stars (10s to 100s M Q ) at very low 
metallicities (Z < 10^ 3 - 75 Z©); a CMB-regulated mode, producing moderate mass (10s of M Q ) stars 
at high mctallicites (Z > 10 _2 5 Z Q at redshift z ~15-20); and a low-mass (a few Mq) mode existing 
between those two metallicities. As the universe ages and the CMB temperature decreases, the range 
of the low mass mode extends to higher metallicities, eventually becoming the only mode of star 
formation. 

Subject headings: stars: formation, cosmology 



1. INTRODUCTION 

Understanding the nature of the first stars in the uni- 
verse is, in principle, a very straightforward problem 
to solve. The initial conditions are de fined by well- 
constr ai ned cosmologica l param eters (e.g.. ISpergel et ail 
(|2007f) ; iKomatsu et all (|2008f) ) and accep ted calcula- 
tions of Big Bang Nucleosynthesis (e.g., ISmith et al.1 
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(|1993t) ). The r elative simplicity of the chemistry o f 
primordial gas (|Abel et all Il997t iGalli fc Pallal 119981) . 
along with powerful numerical methods, have made it 
possible to accurately characterize the formation pro- 
cess of primordial stars from the assembly of their 
dark matter halos through the end of the optically- 
thin regime of the collapsing protostar (|Abel et al.ll2002t 
Bromm et al.ll2002t iBromm fc Loebll2004 lYoshida et al l 
20061: IQ'Shea fc Normanl 120071: IGao et alJ I2007D . ReT 
suits from theory and numerical simulations sug- 
gest that the first stars were 10s to 100s Mp, 
(lAbel et all [2001 lOmukai fc Pallal [200l iTan fc McKeei 
120041: lYoshida et alJ 120061 : IQ'Shea fc Normanl l2007f TA 
direct calculation, however, involves simulating the com- 
plex processes of accretion and radiative feedback, a ca- 
pability currently just beyond the state of the art. In 
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contrast, the challenge of understanding the second gen- 
eration of stars is significantly more complex. The ini- 
tial conditions of the first stars essentially depend only 
on the principles of cold dark matter cosmology and 
the properties of molecular hydrogen, but the initial 
conditions of the second stars require a complete solu- 
tion to the formation and evolution of their predeces- 
sors as well as host of additional chemistry and phys- 
ical processes. The chemical composition and physi- 
cal conditions of second-generation star-forming envi- 
ronments s trongly depend on the exact masses of the 
first stars dHegerfc Wooslevl l2002t iMaeder et "all 120051 : 
iNomoto et al.l l2006t iR ockefeller et alll2006ft and th e me- 
chanics of their supernovae (|Kitavama fc Yoshidal [20051 
iGreif et al.ll200l . 

At the heart of the problem regarding forming the sec- 
ond stars is understanding how the star formation pro- 
cess is altered by the introduction of the first metals in 
the universe, created in the explosions of the first stars. 
The chemical composition of a gas determines the effi- 
ciency with which it can cool radiatively. Cold, metal- 
free gas (T < 10 4 K) cools solely through ro-vibrational 
lines of H2, whose lowest-lying transition has an energy 
equivalent temperature of ~ 512 K, resulting in a min- 
imum achievable temperature of ~ 200 K. The energy 
levels of H2 become thermalized at rather low densities, 
n ~ 10 3 ~ 4 cm~ 3 , where n is the number density, above 
which the cooling rate is only proportional to n, instead 
of n 2 . Numerical simulations have shown that this cre- 
ates a stalling point in the star formation process, where 
the free-fall collapse of metal-free gas comes to a halt as 
the cooling time increases to be s i gnificantly above th e 
dynamical time ([Abel et all 120021 : IBromm et all [2002] ). 
Entry into this 'loitering' phase marks the end of hier- 
archical fragmentation that occurs during the free-fall 
period, a fragmentation that occurs because the tem- 
perature is able to continually decrease with increasing 
density. The final fragmentation mass scale is then set 
by the Jeans mass at the temperature and density cor- 
responding to the point where t he gas can no longe r get 
colder with increasing density ([Larson! I1985L l2005h . In 
the metal-free case, this yields a mass scale of ~ 1000 
Mq, resulting in the high-mass nature of the first stars. 

The addition of metals enhances the cooling rate at 
low temperatures through fine-structure and molecu- 
lar transitions, as well as by continuum emission from 
dust grains. At low to moderate densities (n < 10 9 
cm -3 ), only metals in the gas phase contribute signifi- 
cantly to the coolin g, unless the level of metal enrich- 
ment is very high ( Smith et al.l l20(38h . In the phases 



where H2 cooling is strong (n < 10 cm -3 , an increase 
in the cooling due to metals is first seen only at temper- 
atures less than the 200 K tempera ture minimum tha t 
exists in metal- free star formation ([Smith et al.l l2008h . 
At abundances of Z < 10~ 4 Zq, where metal cool- 
ing is only significant for T < 200 K, the low-density 
evolu tion of the cl o ud is identical to the metal-free 
case ([Omukail 120001 : IBromm et all l200lL lOmukai et all 
12005ft . At higher metallicities, gas-phase metal cool- 
ing becomes strong enough for the collapsing cloud 
to bypass the loite ring phase and to undergo contin- 
ued fragmentation (Bromm et al.l 120011: IBromm fc Loebl 
[2001 ISantoro fc Shulll l2006t iSmith fc Sigurdssonll2"007ir 
If dust grains are present, their influence becomes im- 



portant at very high densities (n > 10 cm ). This 
has been shown to i nduce fragmentation for metallicities 
as lo w as 10~ 5 - 5 Z^ (lOmukai et alj|2005t ISchneider et all 
[2001 [TsuTibe fc Omukail l2006t I Clark et al.ll2008f ). Cal- 
culations bv ISchneider et alJ (|2004l ) predict that up to 
30% of the progenitor mass is converted into dust in a 
pair-instability supernova, but observations of Type II 
supernova in the local universe rem nants have not re- 
turned conclusive evidence of dust ([Green et al.1 120041 : 
iKrause et al.1 12004). In the context of star formation in 
the early universe, the arrival at this critical metallicity, 
Z cr , at which fragmentation occurs beyond the capabili- 
ties of primordial gas, is predicted to be the point where 
the universal mode of star formation shifts permanently 
from the high-mass, solitary mode of the first stars, to 
the low-m ass, multiply-producing mode that is presently 
observed (IBromm fc Loebl 120031: ISantoro fc Shvul 120061 : 
ISmith fc Sigurdssonll2007l) . 

The protostellar collapse of metal-enriched gas 
clouds has been studied extensively with one-zone 
mode ls coupled to lar ge chem ic al networks (Omukail 



ISchneider et alf 120031 : lOmukai et alJ 120051 : 
Schneider et al.l 120061 ). These studies have produced 
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insight into the evolution of the density and temperature 
of collapsing gas clouds with finite metallicities, but 
one-zone models cannot speak to the actual process 
of fragmentation, which requires attention to complex 
cloud geometries that can only be given b y fully 
three -dimensional hydrodynamic simulations ([Larson! 
l2007fl. The first of su ch simulations were carried out 
bv IBromm et al.l (|200l[ ). who included cooling from C, 
N, O, Fe, S, and Si, but not H 2 , finding that clouds 
with metallicities, Z > 10~ 3 Zq are able to fragment 
to mass scales lower than in the metal- free case. These 
simulations had a mass resolution of only 100 M Q and 
were, therefore, unable to investig ate the formation of 
solar-m ass stars. More recently, ISmith fc Sigurdsson! 
(|2007f ) have performed a series of high-resolution simu- 
lations of metal-enriched gas collapse that were able to 
follow the evolution of gas fragments to sub-solar mass 
scales. They, too, find that gas with metallicities, Z 
> 10~ 3 Zq will fragment into multiple clumps, while 
gas with Z < 10~ 4 Zq will produce only one object. 
However, ISmith et al.1 (|2008| ) reported the existence of 
regions in density and temperature that are thermally 
unstable in gas with metallicities as low as 10~ 4 Zq. 
Fragmentation is traditionally thought to happen when 
the cooling time is less than the dynamical time, as 
the gas is able to cool and form perturbations before 
they can be smoothed out by sound waves. When this 
condition, referred to as the fragmentation criterion, is 
satisfied, frag mentation can also be aided by thermal 
instabilities (Field 1965), where a slight decrease in tem- 
perature or increase in density leads to a higher cooling 
rate, causing differences in temperature between pockets 
of gas to grow in a runaway fashi on. In the adaptive 
mesh refinement simulations of Smith fc Sigurdsson 
(2007), grid refinement was performed based on baryon 
and dark matter overdensities and by ensuring that 
the Truelove criterion, lj < Ax, where lj is the local 
Jeans length and Aa; is the grid cell size, was satisfied 
([Truelove et al.l Il997h . However, refinement was not 
performed when the cooling time was less than the 
hydrodynamic time-step. While this was not explicitly 
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TABLE 1 
Simulations 



Run 


IC 


Z (Z Q ) 


Zcol 


Grids 


Cells (xlO 7 ) 


rimax (cm 3 ) 


At co ; (yr) 


rl_rnf 


1 





14.761 


13790 


6.49 


5.96 XlO 11 


- 


rl_Z-6 


1 


10 -6 


14.762 


13015 


6.42 


6.43 xlO 11 


2.9 xlO 4 


rl_Z-5 


1 


10 -5 


14.783 


12889 


6.43 


6.74 xlO 11 


5.7 xlO 5 


rl_Z-4.25 


1 


1Q-4.25 


14.809 


12955 


6.40 


6.75 xlO 11 


1.2 xlO 6 


rl_Z-4 


1 


10 -4 


14.830 


12964 


6.38 


6.22 xlO 11 


1.8 xlO 6 


rl_Z-3.75 


1 


10 -3.75 


14.848 


13003 


6.38 


6.97 xlO 11 


2.2 xlO 6 


rl_Z-3.5 


1 


10 ■ 3 -° 


14.874 


12883 


6.40 


6.89 xlO 11 


2.9 xlO 6 


rl_Z-3.25 


1 


10 


14.936 


12888 


6.36 


3.97 xlO 11 


4.5 xlO 4 


rl_Z-3 


1 


10 3 


15.073 


12581 


6.25 


4.42 xlO 11 


7.9 xlO 6 


rl_Z-2.5 


1 


1Q -2.5 


16.180 


11187 


5.69 


4.70 xlO 11 


3.3 xlO 7 


rl_Z-2 


1 


10 2 


19.481 


8693 


4.57 


1.56 xlO 11 


8.8 xlO 7 


r2_mf 


2 





17.409 


8684 


4.76 


7.97 xlO 11 


- 


r2_Z-4 


2 


10 -4 


17.555 


8509 


4.74 


8.39 xlO 11 


2.5 xlO 6 


r2_Z-3.5 


2 


10 -3.5 


17.654 


8476 


4.71 


8.20 xlO 11 


4.2 xlO 6 


r2_Z-3 


2 


10~ 3 


17.955 


8408 


4.67 


7.30 xlO 11 


9.2 xlO 6 


r2_Z-2.5 


2 


10 -2S 


18.537 


8022 


4.51 


4.41 xlO 11 


1.8 xlO 7 


r2_Z-2 


2 


10~ 2 


20.441 


7194 


4.09 


1.92 xlO 11 


4.4 xlO 7 


r3_mf 


3 





23.885 


7771 


4.25 


1.63 XlO 12 




r3.Z-4 


3 


10" 4 


23.966 


7722 


4.23 


1.65 xlO 12 


6.7 xlO 5 


r3.Z-3.5 


3 


1Q-3-5 


24.122 


7640 


4.21 


1.42 xlO 12 


1.9 xlO 6 


r3_Z-3 


3 


10~ 3 


24.390 


7366 


4.17 


1.29 xlO 12 


4.1 xlO 6 


r3.Z-2.5 


3 


1Q-2-5 


24.732 


7145 


4.13 


1.20 xlO 12 


6.7 xlO 6 


r3.Z-2 


3 


10 -2 


25.028 


7083 


4.09 


6.53 xlO 11 


8.9 xlO 6 


r3_Z-2noCMB 


3 


10~ 2 


25.255 


7424 


4.29 


6.02 xlO 11 


1.1 xlO 7 



Note. — z co i is the redshift at the onset of runaway collapse. The total number of grid 
cells includes those that are covered by child grids at higher levels of refinement. n max is 
the proper maximum baryon number density within the box. At co ; is the time difference to 
runaway collapse from the metal-free case. 



wrong, since the radiative cooling solver iterates with 
time-steps that are no larger than 10% of the cooling 
time, it may have artificially suppressed the growth 
of thermal instabilities that could have formed extra 
fragments. 

We rerun the simulations of ISmith fc Sigurdssonl 
(2007) with an additional refinement criterion which en- 
sures that the hydrodynamic time-step is always less than 
the cooling time on the finest level of resolution. In addi- 
tion, we extend the series of simulations to include lower 
metallicities and to more carefully examine the metallic- 
ity range between 10 -4 Zq and 10 -3 Zq. We also run 
similiar simulations with two extra sets of initial con- 
ditions to confirm the robustness of the results. We de- 
scri be the setup of our simulatio ns and the improvements 
over lSmith fc Sigurdssonl (I2007T ) in §2] In fj3J we calculate 
the gas-phase critical metallicity, evaluate the validity of 
our assumption of optical thinness, and present the re- 
sults of the suite of simulations. In gl we discuss our 
results in the context of star formation at high redshift 
and make the case for an initial mass function (IMF) that 
evolves over cosmic time. We also include a discussion 
of the caveats and limitations of this work. Finally, we 
conclude with a brief summary of the main conclusions 
of this work in fjSJ 

2. SIMULATION SETUP 

We perform a series of 24 primordial star formation 
simulations using the Eulerian adaptiv e mesh refinement 
hydrodynamics + N-body code, Enzo (|Brvan fc Normanl 
119971: lO'Shea et all 120041 ) ■ Excluding our three metal- 
free control runs, the gas in each simulation is homoge- 
neously pre-enriched to some non-zero metallicity. As in 



ISmith fc Sigurdssonl (|2007[ ). we confine the simulations 
to constant metallicities with solar abundance patterns, 
saving the more realistic, and far more complicated, sim- 
ulations of true second-generation star-forming environ- 
ments, with heterogeneous metal-mixing and nonsolar 
abundance patterns, for a future work. 

The nature of the initial c onditions for our simulation s 
are identical to those used in lSmith fc Sigurdssonl (2007). 
The simulation box has a comoving size of 300 h~ 1 kpc, 
with 128 3 grid cells on the top grid and three nested sub- 
grids, each refining by a factor of 2, for an effective top 
grid resolution of 1024 3 cells. The cosmological parame- 
ters have the following values: Qm = 0.3, Sl\ = 0.7, fig 
= 0.04, and Hubble constant, h = 0.7, in units of 100 km 
s _1 Mpc -1 . The po wer spectrum of i nitial density fluc- 
tuations is given by lEisenstein fc Hul (|1999l ). with erg = 
0.9 and n = 1. Refined grids are created during the sim- 
ulations when the baryon (dark matter) density is 4 (8) 
times greater than the mean density at that level. The 
density threshold for refinement decreases at higher lev- 
els. The local Jeans length is resolved by a minimum of 
16 grid cells at all tim es, exceeding the Truelove criterion 
(|Truelove et al.lfl997f ) by a factor of four along each coor- 
dinate axis. In addition, grid refinement occurs whenever 
the cooling time drops below the integration time step of 
the hydrodyna mic solver. This final r efinem ent criterion 
was not used in lSmith fc Sigurdssonl (|2007l ). During the 
simulation, a grid cell is flagged for refinement if one or 
more of any of these criteria are met. 

We perform three sets of simulations. Qualitatively, 
the three sets are the same. They each have the same 
cosmological parameters, box size, and resolution. The 
only difference between them is that their initial condi- 
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tions were created with three unique randomizations of 
the initial density and velocity perturbations. Thus, they 
represent three different realizations of the same prob- 
lem. The first set of i nitial conditions is the one used by 
Smith & Sigurdsson (2007) . The second and third sets 
were initially used bv lO'Shea fc Normanl (|2007l ) and cor- 
respond to the runs named L0_30A and L0_30D in that 
work. For each set of initial conditions, we perform a 
metal-free control run. Excluding the control runs, we 
run 10 simulations using initial condtions Set 1, with 
metallicities ranging from W~ e Z@ to 10 _2 Zq, 5 simu- 
lations using Set 2, and 6 simulations using Set 3, with 
metallicities from \Q~ a Zq to 10~ 2 Z© for Sets 2 and 3. 
The final simulation in Set 3 has a metallicity of 10~ 2 Z Q , 
but unlike the others, excludes the effect of the CMB on 
the cooling of gas. 

We use the second imp lementation o f the o ptically-thin 
metal cooling method of lSmith et alj (|2008[ ). This meth- 
ods uses tabulated cooling functi ons created with th e 
photoionization software, Cloudy (Ferl and et al.l fl998), 
for all elements heavier than He, up to atomic number 
30 (Zn). A solar abundance pattern is used for the metals 
in all of the simulations. Since we only follow the evo- 
lution of the collapsing clouds up to densities of ~ 10 12 
cm~ 3 , we neglect the cooling from dust and assume the 
presence of only gas-phase metals. In §3.3( we calculate 
the optical depth in the collapsing cloud for the most im- 
portant high density, gas-phase metal coolants and find 
that the optically thin assumption holds at all times dur- 
ing simulations. The H/He chem istry is followed ex plic- 
itly during the simu lation as in lAbel et all (|1997l ) and 
lAnninos et al.l (|1997f ). but including reacti ons that ex- 
tend i ts validity to above 10 8 cm -3 . See ISmith et al.l 
(2008) for a full dis cussion of this metho d. We neglect 
cooling from HD, as iBromm et a l. ( 2002) found its con- 
tribution to be negligible in situations where the initial 
ionization is low, su ch as in this w ork. However, it has 
been pointed out by iGloverl (|2008f ) that when the initial 
ionization is high, the HD fraction can become enhanced 
to the point where HD cooling is important. Heating 
from H2 formation is not included in the chemical net- 
work. This may contribute a significant source of heat 
when H2 formation via th ree-body reactions b ecomes im- 
portant at n > 10 9 cm -3 (|Omukai et al .112005?) . However, 
we argue in £13.41 that our results are robust despite this 
shortcoming. The metal cooling data was created with 
the Linux computer cluster Lion-xo which is operated by 
the High Performance Computing Group at The Penn- 
sylvania State University. 

As in the precursor to this work (|Smith fc SigurdssonI 
2007), we do not assume the existence of an ionizing 
UV background. The Population III stars responsible 
for the enrichment of the gas in our simulations have 
died, and we assume any living Pop III stars are most 
likely too distant to contribute significant amounts of 
radiation. A major difference between our simulations 
and those of lSmith fc S igurdssonI (|2007f ) is the choice of 
H 2 cooling rates. We u se the upda ted H 2 cooling rates 
of iGalli fc Palla! ([1998), wh ereas ISmith fc Sigurds"sor] 
d2007D u sed the older r ates of|Lgmi^Siml3l (H984D- The 
rates of [Galli & Palla (199§) are generally lower than 
those of ILepp fc Shul (fl98l . causing the gas to take 
longer to cool in the newer simulations. As a result, the 
moment of runaway-collapse in our simulations is sys- 



tematically delayed from those of ISmith fc SigurdssonI 
(|2007| ). when comparing the simulations in Set 1. 

Each simulation is initialized at z = 99 and runs until 
the point where at least one dense, prestellar core forms 
within a ~ 5x 10 5 Mq dark matter halo, located in 
the center of the simulation box. The simulations are 
stopped when a maximum refinement of 24 levels below 
the top grid has been reached. This corresponds to a 
maximum density of ~ 10 11 cm -3 , or roughly 3 x 10~ 13 
g cm -3 . The simulations were run on DataStar, an IBM 
Power4 machine at the San Diego Supercomputing Cen- 
ter. A summary of the final state of each simulation is 
given in Table [TJ 

3. RESULTS 
3.1. Critical Metallicities 

The gas phase critical metallicity, Z cr , has been es- 
timated analytically by calculating the chemical abun- 
dance required for the cooling time to equal the dynam- 
ical time at the stalling po int for metal-free gas, n ~ 
10 3 " 4 cm' 3 and T ~ 200 K. IBromm fc Loebl (|2003h per- 
formed this e xercise with C and O, bu t excluding cooling 
from H 2 , and iSantoro fc Shulll (1200(1 did so wit h C, O, 
Fe, Si, and including H 2 . ISantoro fc Shulll (|2006f ) consid- 
ered densities above and below the stalling point as well. 
The general consensus from these studies is that Z cr 1=3 
10" 3 - 5 Z Q at n = 10 4 cm' 3 . 

We observe a small amount of variance in the values of 
the density and temperature at which the temperature 
minimum occurs in our three metal-free runs. Since the 
exact value of the critical metallicity depends on the pre- 
cise conditions at the temperature minimum, i.e., tem- 
perature, total density, and H 2 density, we calculate the 
critical metallicity for each of the three runs separately. 
To do this, we define the cooling time in the following 
way: 

nkT 

tcool = 7 ,uh '. Yl Ti (1) 



( 7 -i)[a^ 2 +a; 



etals 



where k is Boltzmann's constant, 7 is 5/3, and A' H2 and 
^'metals are the cooling rates from H 2 and the metals 
in units of [erg s~ 4 cm -3 ]. We use A' to avoid confusion 
with the term, A, which is often expressed in units of [erg 
s~ 4 cm 3 ]. We then set EquationQ]equal to the dynamical 
time, 



tdyn 



3tt 



(2) 



where G is the gravitational constant and p is the mass 
density. We calculate Z cr for the conditions at the tem- 
perature minimum in ea ch metal-free si mulat ion using 
the H 2 cooling rates o f IGalh fc Pallal (1 19981) and the 
metal cooling rates from iSmith et al.l (|2008l ). The results 
are shown in Table [2] The systematic shift toward lower 
values of Z cr compared to previous calculations is most 
likely due to the fact that the minimum temperatures are 
slightly higher than 200 K, where both the H 2 and metal 
cooling rates are higher. In addit ion, cooling from H 2 wa s 
not included in the calculation of lBromm fc Loebl (2003) , 
meaning that more cooling from the metals would have 
been required. There is a correlation between Z cr and the 
collapse redshift of the simulation, with the highest value 
of Z cr coming from the highest collapse redshift, despite 
the lack of such a correlation for n, nn 2 , or T. However, 
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TABLE 2 
Critical Metallicities 



Run 


n [cm 3 ] 


n H2 [cm 3 ] 


T 


log(Z cr /Z ) 


rl_mf 


6.89X10 3 


3.42 


283 


-4.08 


r2_mf 


3.64X10 3 


1.86 


214 


-3.90 


r3_mf 


1.19X10 4 


6.53 


260 


-3.85 



Note. — Critical metallicities calculated for each 
mctal-frcc simulation, n, nu 2 > an d T are the proper 
number density, H2 number density, and temperature 
at the temperature minimum where H2 becomes ther- 
malized in each of the metal-free simulations. The final 
column is the log of the metallicity required to equate 
the cooling time to the dynamical time for the condi- 
tions listed. 

there is a correlation between the H2 fraction and col- 
lapse redshift, which has been observe d in the Pop III 
simulations of lO'Shea fc Normanl (|2007l) . They find that 
the higher H2 fractions result from the generally warmer 
gas in halos that collapse at higher redshifts, which is 
simply a function of the linear dependence upon redshift 
of the virial temperature. This may be the dominant 
factor, but with only 3 data points it is unclear whether 
the observed trend is even significant. In addition, the 
artificial nature of our initial conditions may make this 
finding inapplicable to the real world. 

We see a trend with metallicity and collapse 
redshift that is si m ilar t o what was reported by 
iSmith fc Sigurdssonl ([2007), where simulations with 
higher metallicities reach the runaway collapse phase ear- 
lier. We define At co i as the difference in time to runaway 
collapse between a simulation with non-zero metallicity 
and the metal-free run with the same initial conditions. 
An increase in metallicity by 0.5 dex results in an increase 
in At co i by a factor of approximately 1.3 to 4. Similarly, 
when compared within a set, the simulations collapsing 
later have a higher number of total grids and grid cells 
in their final output, since the lower-density envelope gas 
has had more time to evolve and reach higher densities. 

3.2. Radial Profiles 

In Figure [TJ we show projections of mass- weighted 
mean number density for the central 0.5 pc surrounding 
the point of maximum baryon density for the final out- 
put of all runs in Set 1. For the runs with metallicities 
near or below Z cr (lO -4 ' 08 .^ for this set of simulations) , 
the central cores appear quite round and show no clear 
signs of forming more than one object. In the metallicity 
range from 10~ 3 ' 75 Z Q to 10~ 3 25 Z Q , the cores appear in- 
creasing asymmetric, with at least one additional density 
maximum present. However, at metallicities at or above 
10~ 3 Zq, the cores return to a more spherical shape with 
only a single density maximum. 

In Figure [21 we plot spherically averaged, mass- 
weighted values of the density, temperature, enclosed 
mass, and radial velocity as a function of radius from 
the point of maximum density for all simulations in Set 
1. On large scales, the density profiles follow an r~ 2 - 2 
power law. To highlight the difference in density be- 
tween each run, we plot in panel A of Figure [2] the value 
of (n x r 2 ), instead of simply n. We choose to scale the 
density by r 2 instead of r 2 2 because it is easier to glean 
the true density from the figure. Over most of the plot- 



ted range, the run with the lowest metallicity has the 
highes t density. In the isothermal collapse model of lShul 
(|1977l ). the accretion rate is proportional to the cube of 
the sound speed, or T 3 / 2 . Figure [2)3 shows that while 
isothermality does not really apply, there is a clear cor- 
relation between the temperature and density. The runs 
with the highest metallicity, and subsequently the coldest 
gas, are the least dense. In addition, within individual 
runs, an increase in the gas temperature is matched by 
an increase in the density. 

The instantaneous accretion rate at a given radius is a 
function of the infall velocity at that position. Although 
there are some exceptions, the correlation between tem- 
perature/sound speed and infall velocity, with higher 
temperatures/sound speeds corresponding to higher ve- 
loc ities, generally holds. Th is was also found to be true 
bv lO'Shea fe Normanl {2007) in their simulations of Pop 
III star formation. We find that the inflow is roughly 
transonic throughout the entire density range, showing 
that even though the collapsing clouds are not isothermal 
spheres, their accretion rates are still largely regulated by 
the sound speed. 

3.3. Optical Depth 

The fragmentation of collapsing gas depends sensi- 
tively on its thermal evolution, which is effectively a mea- 
sure of the radiative cooling properties of the gas as a 
function of density and temperature. Since our radiative 
cooling method, described above, assumes optical thin- 
ness, the results of our simulations are only valid where 
r <C 1. The optically thin assump tion begins to brea k 
down at densities n > 10 10 cm -3 (jOmukai et alJl2005Q . 
The optical depth at a frequency v is expressed as 

t v = J k v p d£, (3) 

where k„ is the opacity, p is the mass density, and d£ is 
the distance traveled by a photon. At a given density, the 
cloud has a characteristic size r(p), shown in Figure [2JA.. 
If a photon emitted from a region with density p is able 
to travel a distance ~ r(p) without being re-absorbed, 
then it will not alter the thermal evolution of the cloud, 
and the assumption of optical thinness holds. Using the 
Cloudy software, we calculate values of the absorption 
coefficient, a v = n v p, as a function of density, metallicity, 
temperature, and frequency. Equation [3] then takes the 
form 

t v =a v {p, Z,T)r(p). (4) 

For metallicities Z < 10~ 2 Z Q and in the absence of dust, 
metal cooling is dominated by fine-structure transition s 
of Oi and Fen (|Santoro fc Shulll[200l ISnnth et al.l[200l . 
In Figure [3l we plot the optical depth from Equation [4] 
for the Fen line at 25.99 pm and the Oi line at 63.18 pm 
as a function of density using the spherically-averaged 
densities and temperatures shown in Figure [2j At low 
energies and in the absence of dust grains, the largest 
contributor to the opacity is free-free absorption, which 
is dominated by H for low metallicities. As a result, the 
opacity in the energy range of interest is essentially in- 
dependent of metallicity for metallicities less than solar. 
When the metals become a considerable fraction of the 
total gas content, above Z ~ 1OZ or so, the opacity at 
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Fig. 2. — Spherically-averaged, mass-weighted quantities as a function of radius from the point of maximum density for 7 of the 10 
simulations in Set 1. A: number density normalized to an r~ 2 power law; B: temperature; C: enclosed gas mass; D: radial velocity. In 
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Fig. 4. — The number of bound clumps found within a sphere of 
radius 5 pc, centered on the point of maximum density, as a func- 
tion of metallicity for all of the simulations performed. The grey 
bar in the bottom panel denotes the run where the temperature 
floor created by the CMB was not included. In each panel, the 
vertical dashed line represents the estimated value of Z cr for each 
set from Table |2] 

low energy for dust-free gas begins to grow with metal- 
licity. In Figure El the apparent decrease in the optical 
depth at the highest metallicities is due to the flattening 
of the density profile at small radii. As such, the final 
point for each curve in Figure [3] should be ignored. For 
values inside the point where the density profile flattens, 
it would be more reasonable to use the radius at which 
the density turns over for the calculation of r. Fortu- 
nately, r remains significantly less than 1 throughout all 
of our simulations, peaking at roughly 0.25 at the high- 
est densities for the runs with the lowest metallicities. 
The lower optical depth in the higher metallicity runs 
is due to the smaller characteristic size of the core, as 
seen in Figure [2JA It should also be noted that some 
molecules are missing from the Cloudy dust-free chem- 
istry network that could contribute to the opacity, such 
as TiO. As such, the optical depths calculated could be 
slightly higher for the high metallicity simulations. How- 
ever, we show in £)3.4I that fragmentation in our simula- 
tions occurs well before this becomes a concern. 

3.4. Fragmentation 

In order to quantify the degree of fragmentation within 
each run, we employ a clump finding algorithm to search 
for b ound clumps within 5 pc of the density maximum. 
As in lTruelove et all (|1998f ). we define a clump, there re- 
ferred to as a fragment, as "the mass contained between a 
local density maximum and the lowest isodensity surface 
surrounding only that maximum," to quote that work. 
We begin by identifying all grid cells within a sphere of 
radius 5 pc, centered on the point of maximum density. 
We then create density contours using all the cells within 
the sphere. In density space, the first contour spans the 
entire range of density within the sphere, effectively cre- 
ating one large contour. The contour becomes the par- 
ent clump of all other clumps that will be found as the 
process continues. On the second iteration, we create 
contours with the same maximum as before, but with 
the minimum increased by 1/4 dex. If more than one 
contour exists, those groups of cells become child clumps 
of the group made by the previous iteration. The pro- 
cess continues in a recursive fashion, creating groups of 
cells based on contours of increasing minimum density, 



TABLE 3 

Fragmentation Properties 



Run 


Clump Mass [M ] 


n frag [cm 3 ] 


rl Z-3 5 


KA A A QQ 




rl_Z-3.25 


22.1 


1.55 xlO 6 




5.46, 1.40 


8.73 xlO 6 


r2_Z-3 


6.26, 1.89 


2.75 xlO 9 


r2_Z-2.5 


95.6, 27.7 


2.76 xlO 4 


r3.Z-3.5 


54.4, 4.98 


1.55 xlO 5 


r3.Z-2 


289, 105 


1.55 xlO 4 


r3.Z-2noCMB 


11.7 


2.76 xlO 4 




4.38, 4.06 


4.91 xlO 5 




32.6, 0.49 


8.73 xlO 5 



Note. — Number density at which fragmenta- 
tion occurs within simulations where multiple bound 
clumps are identified and the masses of the clumps 
associated with that episode of fragmentation. The 
fragmentation density is taken to be the minimum 
density within a clump, taken from the earliest data 
output where that clump is identified. 

identified only within the cells of the parent clump. The 
clump finding ends when minimum countour density has 
reached the constant maximum. Effectively, we create a 
family tree of clumps, with the very first group as the 
trunk of the tree. During the clump finding process, 
a child clump is only kept if it is gravitationally self- 
bound, or has children of its own that are bound. When 
determining whether a clump is bound, we consider the 
thermal energy of the gas as well its kinetic energy with 
respect to the bulk center of motion. 

In Figure 01 we plot a histogram of all the number of 
clumps found in each run within the 5 pc sphere, as a 
function of the metallicity of the run. For this plot, we 
only include bound child clumps with no children of their 
own. When enlarging the radius of the sphere from 5 pc 
to 10 pc, no additional clumps were found in any of the 
runs. Figure [4] confirms what is seen in Figure [T] In all 
runs with metallicities below Z cri only a single bound 
clump is found. As the metallicity increases, the number 
of clumps increases, then decreases back to only a single 
clump for the highest metallicities, with the exception of 
r3_Z-2, which has 2 bound clumps. The range of metal- 
licities where fragmentation occurs is consistent between 
Sets 1 and 3, but offset by 0.5 dex toward higher metallic- 
ities for Set 2. It is not clear what causes this offset, but 
the qualitative trend of increasing and then decreasing 
number of clumps exists in all 3 sets. It is also worth- 
while to note that runs rl_Z-4, rl_Z-3.75, and r2_Z-3.5, 
while slightly above Z cr , do not show fragmentation. 

During the simulations, the Enzo code creates a snap- 
shot of the entire box each time the maximum level of 
refinement increases. This provides us with multiple data 
outputs as the central density increases during runaway 
collapse. We ran the clump finder on all data outputs 
created during this period, searching for the first data 
output in which multiple clumps arc found. Within this 
output, we take the minimum density within the clump 
to be the density at which fragmentation occurred. In 
Table O we list the fragmentation densities for all sim- 
ulations in which multiple clumps were found. In runs 
with 2 clumps, the fragmentation density is the same 
for both clumps, since it represents the lowest density 
contour for which the two objects are separate. In run 
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rl_Z-3.25, which has 3 clumps, the core initially frag- 
mented into 2 clumps at n = 1.55 x 10 6 cm' 3 . One data 
output later, one of those clumps fragmented again at 
n = 8.73 x 10 6 cm -3 . A similar thing occurred in run 
r3_Z-2noCMB with 4 clumps forming initially and one of 
those fragmenting again later. 

At densities of roughly 10 9 cm -3 , H2 formation via 
three-bo dy reactions begin to rapidly increas e the H 2 
fraction (|Abel et all 120021 : iBromm et all 12002ft . In our 
simulations, we find the H2 fraction to be roughly con- 
stant at ~ 10' 3 until the density reaches ~ 10 9 cm' 3 . 
At densities of a few x 10 10 cm -3 , the H2 fraction has 
passed 10% and continues to rise. The onset of rapid 
H2 formation occurs slightly earlier (within a factor of a 
few in density) in the runs with higher metallicity due to 
the T' 1 dependence of the three-body rate coefficients. 
Our chemical network lacks heating from H2 formation 
(~ 4.4 eV per reaction), which may significantly raise 
the temperature of the gas when three -body reactions 
become important (jOmukai et al.l 12005ft . We test this 
by calculating the ratio of the heating rate produced by 
H2 formation to the total rate of cooling (without H2 
formation heating) for each of the simulations. We use 
the rate coefficients k^ (H + H + H — > H2 + H) an d ke 
(H + H + H 2 -> H 2 + H 2 ) from lPalla etaD (fl983h and 
assume an energy injection of 4.4 eV per reaction. We 
also include a cooling term representing the reverse of 
the above two rea ctions, using coefficients k$ and k^ of 
iPalla et all (|1983ft . For metallicities Z < 1O' 3 - 5 Z , the 
ratio of the H2 formation heating to the total cooling is 
0.1 at n ~ 10 9 cm -3 and 1 at n ~ 5 x 10 9 cm' 3 . For 
Z = 10' 3 Z o , this ratio is 0.1 at n — 9 x 10 7 cm' 3 and 
1 at n ~ 6 x 10 s cm' 3 . For Z > 10' 25 Z Q , this ratio is 
0.1 at n ~ 3 x 10 7 cm' 3 and 1 at n ~ 10 8 cm' 3 . For run 
r3_Z-2noCMB, which has the CMB removed, this ratio 
is 0.1 at n ~ 10 7 cm' 3 and 1 at n ~ 10 8 cm' 3 . For 
simulations with equivalent metallicities, the above ratio 
varies maximally by a factor of a few due to minor differ- 
ences in the temperature at a given density. In all cases 
but one, fragmentation occurs well before H 2 formation 
heating becomes important. The lone exception is run 
r3_Z-3, where fragmentation occurs at n ~ 3 x 10 9 cm' 3 
and H2 formation heating matches the total cooling at 
n ~ 10 9 cm' 3 . However, the core in this simulation does 
not fragment until it has already begun to reheat due to 
adiabatic compression, so it is not clear what effect the 
inclusion of this missing heating term would have had. 

In Figures [5] and [6l we plot the number density vs. gas 
temperature for the final output of each simulation. Due 
to the self-similar nature of the collapse, Figures [5] and [6] 
can also be used to understand the evolution of the cen- 
tral core throughout the collapse. In all the runs with 
Z < Z cr , the cooling is too low to prevent the temper- 
ature from rising at the H2 thermalization density, n ~ 
10 4 cm -3 . Therefore, the minimum fragmentation mass 
for these runs, set by the Jeans mass as the temperature 
minimum, is well over 1,000 Mq, which is nearly equiva- 
lent to the total enclosed mass. Even though runs rl_Z-4 
and rl_Z-3.75 are above Z cr , the additional cooling pro- 
vided by the metals is not sufficient to significantly lower 
the minimum fragmentation mass. For runs rl_Z-3.5 and 
rl_Z-3.25, the more efficient cooling lowers the minimum 
fragmentation mass to just over 100 Mq, which is ap- 
proximately a factor of a few lower than the total mass 



within 1 pc. 

For the runs with the highest metallicities, as in runs 
rl_Z-2.5 and rl_Z-2, the gas cools all the way to the tem- 
perature of the CMB. The cooling proceeds so efficiently 
that the gas has not had sufficient time to reach high den- 
sities before hitting the temperature floor of the CMB. 
Fragmentation can only continue as long as the tem- 
perat ure decreases with increasing density (|Larsonlll985l 
2005). Although the temperature decreases slightly in 
runs rl_Z-2.5 and rl_Z-2 for densities greater than 10 3 
cm' 3 , the temperature minimum is effectively at n = 
10 3 cm -3 , where the gas reaches the CMB temperature. 
Near the CMB temperature, the value of the cooling rate, 
A, effectively becomes (A(T) - A(T C mb))- Therefore, 
when the gas reaches the CMB temperature, the cooling 
rate drops to zero, and the cooling time becomes infi- 
nite. The gas cloud becomes extremely thermally stable, 
preventing further fragmentation. 

To verify that the CMB is indeed suppressing frag- 
mentation, we run one simulation, r3_Z-2_noCMB, with 
the CMB temperature floor removed. We choose ini- 
tial conditions Set 3 for this exercise since it has the 
highest CMB temperature at the redshift of collapse 
and should therefore show the greatest contrast with 
the CMB removed. In Figure [7l we show mass- weighted 
mean number density projections of the central 5 pc for 
runs, r3_mf, r3_Z-2, and r3_Z-2_noCMB. Run r3_Z-2 has a 
much clumpier structure than its metal- free counterpart, 
even with the CMB temperature floor present. How- 
ever, when the temperature floor is removed, the gas is 
able to collapse into a long, thin filament with far more 
small-scale structure. As shown in Figure |4l we find the 
most bound clumps in this run (4 within 1 pc of the 
density peak and 1 more within 5 pc) . In Figure [8l we 
show mass-weighted mean temperature projections for 
the same runs as Figure [7] We overlay contours of pro- 
jected mean number density of 10 cm' 3 . Figure [6] shows 
that this is the approximate density at which the cloud 
in run r3_Z-2 first reaches the CMB temperature. In 
run r3_Z-2noCMB, it took -400,000 years for the central 
density to increase from 10 4 cm' 3 to 10 5 cm' 3 , which 
is similar to the timescale in a dynamical collapse. In 
run r3_Z-2, the equivalent change in density took —1.9 
million years, indicating that cooling to the CMB tem- 
perature has indeed ended free-fall collapse. 

In the projections of run r3_Z-2 in Figure [H the 2 
largest contours roughly represent the 2 bound clumps 
found within the 5 pc radius sphere. The gas within the 
clumps has a very uniform temperature, as its cooling 
has been abruptly halted at the CMB temperature (—71 
K). The 2 clumps have masses of roughly 100 and 300 
Mq. It is unlikely that they will fragment further, since 
the Jeans mass of each clump is nearly equivalent to its 
total mass. At a temperature of 71 K, sound waves will 
travel — 1.7 pc in the 1.9 Myr required for any additional 
fragments to increase in density by an order of magnitude 
once they have cooled to the CMB temperature. Any ad- 
ditional fragments that might possibly condense out of 
lower densitiy gas must be approximately this large in 
order to collapse, suggesting they will also be very mas- 
sive. It is interesting to note that the two bound clumps 
in the simulation are also about this size. 

In contrast, the gas inside the contours of run r3_Z- 
2noCMB shows significantly more structure in tempera- 
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Fig. 5. — Mass-weighted, average temperature as a function of number density for all runs in Set 1. The colors are the same as in Figure 
[2] including the runs with metallicities Z = 10 _4 25 Z Q (dashed-yellow), 1O _3 - 75 Z (dashed-green) , and 1O _3 ' 25 Z (dashed-blue) . The 
thin, black, dashed lines indicate lines of constant Jeans mass in Mq. The horizontal, blue, dashed line denotes the temperature of the 
CMB at z = 19, the approximate rcdshift of collapse for runs rl_Z-2.5 and rl_Z-2. The central cores in these two runs were both able to 
cool to the temperature of the CMB. 
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Fig. 6. — Mass-weighted, average temperature as a function of number density for all the runs in Set 2 (left) and Set 3 (right). For 
both panels, the metallicities are Z = (solid- black), 1CT 4 Z© (green), lO~ 3 - 5 Z e (blue), 1CT 3 Z© (purple), W~ 2 - 5 Zq (dashed-black) , and 
10 — 2 Z© (dashed-red). In the bottom panel, the dashed-orange line corresponds to run r3.Z-2_noCMB, with Z = 10~ 2 Z e , but with the 
CMB temperature floor removed. All curves in Figures [5] and [6] with the same colors refer to simulations with the same metallicities. The 
thin, black, dashed lines indicate lines of constant Jeans mass in Mq. The horizontal, blue, dashed lines denote the temperature of the 
CMB at z = 19 (left) and 25 (right). 
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ture. The cold knots in run r3_Z-2noCMB seen in Figure 
M correspond to the high density regions seen in Figure El 
4 of the 5 bound clumps are Jeans unstable, with masses 
of 4, 4, 12, and 33 M Q . The fifth, with M ~ O.5M , 
is approximately 1/6 of its Jeans mass. The most mas- 
sive clump is about 5 times more massive than its Jeans 
mass. 

4. DISCUSSION 

We have shown that fragmentation occurs within a col- 
lapsing cloud when the metallicity is above the critical 
metallicity. The exact value of the critical metallicity 
required to prevent an increase in temperature at the 
stalling point of H2 varies slightly from halo to halo. 
Within our three sets of initial conditions, the values of 
Z cr are correlated to the collapse redshift of the metal- 
free runs, with the highest Z cr corresponding to the high- 
est redshift. However, it is unclear whether this is sig- 
nificant. If the metallicity is only marginally higher than 
Z cr , fragmentation is unlikely to occur, as the increase 
in temperature at the H2 stalling point is only delayed 
momentarily. Thus, it is unlikely that a sharp transition 
in star formation mode occurs at just the moment when 
the critical metallicity is reached. 

Fragmentation is suppressed when the metallicity is 
high enough such that the gas is able to cool to the tem- 
perature of the CMB when the central density is still 
relatively low. We confirm that the CMB is responsible 
for the observed suppression of fragmentation by running 
an identical simulation without the CMB. In this simu- 
lation where the CMB is absent, we find more bound 
clumps than in any other of the runs in this study. 2 
bound clumps were found in run r3_Z-2, where the gas 
was able to cool to the CMB temperature. However, both 
of these clumps were quite massive (M > IOOMq), and 
we showed in H3A\ that is is unlikely that they will frag- 
ment into smaller objects. We observe a small amount 
of variance in the metallicity range in which fragmen- 
tation occurs that does not appear to be related to the 
CMB. Just as the exact value of Z cr seems to vary from 
halo to halo, we suspect that the range of metallicities 
where fragmentation occurs will also be influenced by the 
individual properties of a halo as well as its particular 
evolution. 

The mass scale of collapsing clumps can be estimated 
from the Jeans mass at the end of the cooling phase. 
This implies the existence of three distinct metallicity 
regimes for star formation. In the first regime, Z < Z cr , 
which we refer to as the 'primordial' mode, metals do 
not provide enough additional cooling to allow the gas 
temperature to continue to decrease monotonically with 
increasing density when the core reaches the H 2 thermal- 
ization density. In this case, the collapse proceeds in a 
similar way to the metal-free scenario, resulting in the 
formation of a single, massive object. There is also the 
potential for stars forming in this mode to be somewhat 
less massive than the very first stars. The fragmentation 
mass scale for extremely low metallicity gas {Z <C Z cr ) 
may be lowered though compression by shocks from Pop 
III supe rnovae and enhanced cooling from HP in relic H II 
regions (jMackev et al.ll2003t I Johnson fc Bromml l2006). 

At the other extreme, we define Zcmb as the metallic- 
ity at which the gas can cool to the CMB temperature. 
When Z ^> Zcmb, the cloud-core will efficiently cool to 



the temperature of the CMB when the central density is 
still relatively low. In this scenario, fragmentation is lim- 
ited by cooling rapidly to the CMB temperature, as the 
mass scale is determined by the Jeans mass at the den- 
sity when the core first reaches the CMB temperature. 
We refer to this as the CMB-regul ated star formation 
mode, similar to iTumlinsonl (|2007bf h As fragmentation 
is severely limited in this mode, these stars will most 
likely be more massive on average than the characteris- 
tic mass of stars forming today. 

Finally, our simulations have shown that there exists a 
special range in metallicity, Z cr < Z < Zcm b , where the 
core does not reheat at the metal- free stalling point, but 
also cannot cool all the way to the CMB temperature. 
The minimum temperature is set only by the balance of 
radiative cooling and adiabatic heating. The mass scale 
is not regulated externally by the CMB, but rather inter- 
nally by the metallicity-dependent gas-cooling. Hence, 
we term this the metallicity-regulated star formation 
mode. This mode produces the lowest mass stars of the 
three modes mentioned. 

The CMB-regulated star formation mode creates a 
means by which a higher number of massive stars are 
formed in very early universe, when the CMB tempera- 
ture was much higher. As the universe evolves, the CMB 
temperature will slowly decrease, which will increase the 
metallicity required to reach the CMB temperature, re- 
ferred to here as Zcmb- The decrease in the CMB tem- 
perature also means that the fragmentation mass scale 
will be lower at the point where the gas reaches the tem- 
perature floor. Thus, the characteristic mass of stars 
produced by the CMB-regulated mode will slowly de- 
crease with time. This behavior is in agreement with 
the mo del of an IMF that evolves with redshift formu- 
lated bv lLarsonl (|1998( ). As the metallicity threshhold for 
the CMB-regulated mode advances to higher metallicity, 
the range of operation of the metallicity-regulated mode 
will extend to take its place. We lack sufficient data in 
this study to predict the evolution of Zcmb with red- 
shift. However, in a paper to follow, we will map out the 
evolution of Zcmb vs. z with additional simulations col- 
lapsing as much lower redshifts. Observations of nearby 
star-forming clouds show that the minimum achievable 
tempe r ature in the local universe is roughly 10 K (e.g., 
lEvansI (|1999h ). This implied that the CMB-regulated 
star formation mode is in operation up to 2 ~ 2.7, 
at the absolute latest. A growing amount of evidence 
has been pres ented that t h e stel la r IMF evo l ves wit h 
redshift (e.g iFardal et all (I2007T): ITumlinsonl (l 2007bf) : 
Ivan Dokkuml (|2008MDavlT2 008): Wilkins et al. (2008^. 
Interestingly Ivan Dokkum (|2008l ). iDavel (|2008[ L and 
Wilki ns~et alj (|2008h report evidence from high redshift 
that the IMF may deviate from the standard Salpeter 
IMF, favoring higher mass stars for z > 2 - 4. 

If dust is present in the very early universe, this 
would extend the range of the metallicity-regulated 
star forma tion mode to metal li cities as low as Z ~ 
lQ- 5 - 5 ^ (fOmukai et all 120051: ISchneider et all 120061 : 
iTsuribe fc Omukail 120061 : IClark et all l2008fl . The exis- 
tence of low mass, hyp er metal poor stars, H E0107-5240 
(IChristlieb et alJl2002f ). and HE1327-2326 (jFrebel et all 
2005), both with [Fe/H] < -5, may provide evidence of 
this. Both of these stars show extremely enhanced C and 
O abundances, which would make their effective metallic- 
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Fig. 7. — Projections of mass-weighted mean number density along the x (left), y (center), and z (right) axes for the final output of runs 
r3.mf with zero-metallicity (top), r3.Z-2 with Z = 10~ 2 Z Q (middle), and r3.-2noCMB with Z = 10~ 2 Z Q and the CMB temperature floor 
removed (bottom). Each projection is centered on th e location o f maximum density in the simulation box and has a width of 5 pc proper. 
The images were made with the YT analysis toolkit (Turk 2008, yt.enzotools.org). 



ities (in terms of the radiative cooling ability of gas with 
that a bundance) much hi gher. However, it is pointed 
out by ITumlinsonl ( 2007a. b), that the abundance pat- 
terns of these stars are best recreated by a scenario in 
which the C and O enhancement comes via binary mass 
transfer from an intermediate mass AGB star, meaning 
the stars are truly metal poor. In that case, these two 
low-mass stars would like ly require dust in o rder to form 
at such low metallicity. ITumlinsonl (|2007bl ) also shows 
that such a high fraction of carbon enhanced metal poor 
stars (CEMPs) requires a higher than normal fraction 
of more massive st ars that go through the AGB phase. 
ITumlinsonl (|2007bO claims that the evolution of CEMP 
fraction with metallicity (with a higher CEMP fraction 
at lower metallicity) already shows the influence of the 
CMB on the IMF. If a CEMP star and its binary compan- 



ion were formed from gas at the same metallicity, there 
would have to be process at work that would prevent the 
dust cooling fragmentation from fo rming only low mass 
stars. The star formation models of lOmukai et al.l (2005) 
do not indicate that the dust cooling phase that induces 
low mass fragmentation is able to reach the CMB tem- 
perature for Z < 10 Zq. This may simply imply that 
metal mixing from the first supernovae is higly heteroge- 
neous, allowing stars to form simultaneously with largely 
different abundances. 

In this work, we study only gas clouds with solar abun- 
dance patterns. Most likely, the first metals in the uni- 
verse will not have solar abundance patterns. However, 
from the perspective of simulating metal-enriched star 
formation, it is not the specific elemental abundances of 
a gas cloud that are important, but rather the total cool- 
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roper. The solid lines show contours of projected mean number density of 10 4 cm -3 . The images were made with the YT analysis toolkit 
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ing rate produced by the gas. Therefore, given that it is 
the sum of the metals that is important, and not the 
abundance pattern, the results of this work are robust in 
spite of the fact that the abundance patterns used are 
likely to be incorrect. 

There appears to be some discr e pancy between our 
results and those of Uappsen et alJ (|2007f ) . who see no 
evidence of fragmentation induced by gas-phase metal 
cooling. This is potentially resolved by the fact that in 
that work, the gas collapse is only strictly followed up 
to densities of 5xl0 2 cm" 3 before sink particles are cre- 
ated. Figures \5\ and [6] show very little difference in the 
thermal structure of the gas for densities below 5xl0 2 
cm" 3 . Additionally, our clump finding algorithm found 
only a single bound clump within every simulation when 



the analysis was performed on data outputs that were 
made when the maximum density was only ~ 10 3 cm" 3 . 
Multiple clumps were only found when the clouds had 
reached somewhat higher densities. Finally, the simula- 
tions in this work began with cold, neutral gas, whereas 
their simulations began with hot, ionized gas. Had their 
simulations been run to higher densities, any dissimilari- 
ties might also be due to using different initial conditions. 

5. CONCLUSION 

We have performed a scries of high resolution simula- 
tions of metal-enriched star formation using cosmologi- 
cal, Pop III style initial conditions, and assuming fully 
homogeneous metal enrichment. We have shown that our 
results apply to more than a single star forming region 
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by using 3 different sets of initial conditions with iden- 
tical cosmological parameters and resolution, but with 
3 unique random seeds with which to create the initial 
perturbations in the density and velocity fields. From 
the results of these simulations, the main conclusions of 
this work are: 

1. Fragmentation does not occur when the metallic- 
ity is only slightly above Z cr , since this only leads to a 
small delay in the onset of the loitering phase that is 
brought on by a decrease in the efficiency of H2 cooling. 
The density at which the temperature begins to increase 
with increasing density is only marginally higher than in 
the metal-free case, and therefore does not lead to a sig- 
nicant lowering of the minimum Jeans mass. Within our 
simulations, Z cr is roughly 10~ 3 9 Zq. We find that frag- 
mentation does not occur until the metallicity is roughly 
0.5 dex above Z cr . 

2. Fragmentation is suppressed when the metallicity 
is high enough such that the gas is able to cool to the 
temperature of the CMB when the density of the col- 
lapsing cloud is still relatively low (n ~ 10 4 cm -3 ). The 
Jeans mass at the density and temperature at which the 
cloud first reaches the CMB temperature sets the mini- 
mum fragmentation mass within the cloud, and as such 
only massive clumps are able to form. 

3. Metal-enriched star formation occurs in three modes 
that are separated by two metallicity thresholds, Z cr and 
Zcmb- Z cr is the conventional critical metallicity and 
Zcmb is the metallicity at which the gas is able to cool to 
the temperature of the CMB at a given redshift. At the 
approximate collapse redshifts of our simulations, z ~ 20, 
Zcmb is between 10 -3 and 10 _2 5 Z Q . The three modes 
of star formation are: 

• The primordial mode (Z < Z cr ) - The additional 
cooling provided by the metals is not enough to sig- 
nificantly alter the thermal structure of the cloud 
relative to the metal-free case. No fragmentation 
occurs and the star will have a mass similar to a 
Pop III star. 

• The metallicity-regulated mode (Z cr < Z < 
Zcmb) - Metal cooling is high enough to allow the 
cloud to continue to cool past the metal-free loi- 
tering phase, but not high enough to cool it to the 
temperature of the CMB. The minimum fragmen- 
tation mass is set at lower temperatures and higher 
densities than in the primordial case and fragmen- 
tation into multiple objects occurs. Based on the 
masses of clumps formed within our simulations, 
stars forming in this mode could have masses of 



the order of a few Mq or less. 

• The CMB-regulated mode (Z > Z C mb) - Frag- 
mentation is suppressed when the cloud is able 
to cool to the CMB temperature, as is described 
in point 2 of the conclusions. At minimum, 
stars forming in the CMB-regulated mode will be 
more massive than stars forming in the metallicity- 
regulated mode. Zcmb will increase as the CMB 
temperature lowers with time. As such, the masses 
of stars forming in the CMB-regulated mode will 
slowly decrease. This mode will vanish altogether 
when the CMB temperature reaches the observed 
minimum temperature of nearby molecular clouds 
(T ~ 10 K at z ~ 2.7). 

4. By demonstrating that the CMB can suppress 
fragmentation, we have provided a key co nceptual piece 
CMB- I MF hypoth e sis (iTumlinsonl l2007bD . As pointed 
out by ITumlinsonl (|2007bl ). an IMF that evolves with 
redshift, producing more massive stars in the past, has 
consequences that may already be testable by observa- 
tions. 
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